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The large acceptance and high momentum resolution as well as the significant particle 
identification capabilities of the NA49 experiment |]| at the CERN SPS allow for a broad 
study of fluctuations and correlations in hadronic interactions. 

In the first part recent results on event-by-event charge and (pt) fluctuations are pre- 
sented. Charge fluctuations in central Pb+Pb reactions are investigated at three different 
beam energies (40, 80, and 158 AGeV), while for the (p t ) fluctuations the focus is put on 
the system size dependence at 158 AGeV. 

In the second part recent results on Bose Einstein correlations of h~h~ pairs in mini- 
mum bias Pb+Pb reactions at 40 and 158 AGeV, as well as of K + K + and K~K~ pairs 
in central Pb+Pb collisions at 158 AGeV are shown. Additionally, other types of two 
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particle correlations, namely 7rp, Ap, and AA correlations, have been measured by the 
NA49 experiment. Finally, results on the energy and system size dependence of deuteron 
coalescence are discussed. 

1. Fluctuations 

1.1. Charge fluctuations 
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Figure 1. $ q as a function of the number of 
charged particles in the acceptance window 
Ay for central Pb+Pb collisions. 
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Figure 2. A$ q as a function of the num- 
ber of charged particles in the acceptance 
window Ay for central Pb+Pb collisions. 



It has been suggested that charge fluctuations might be sensitive to the presence of a 
quark qluon plasma phase @||. A suitable observable, however, must take into account 
the effect of impact parameter fluctuations, acceptance effects, the total net charge of the 
reaction system, and charge conservation. For this purpose a generalized <£> x measure [Q 
can be employed: 

$ x = y/(Z 2 }/(N) - \IJ 2 with:Z = EiIx^ and Zi = Xi -x (1) 

Here N is the event multiplicity, while x is the quantity to be studied. The overline denotes 
the average over a single particle inclusive distribution, whereas (■ ■ •) means the average 
over all events. In the case of charge fluctuations, x is chosen as the electric charge q of a 
particle ||. $ q can vary between two extreme cases: For independent particle emission $ q 
is equal to 0, while local charge conservation would imply $ q = —1. Figure |1] shows the 
measured $ q values in central Pb+Pb reactions^] at 40, 80, and 158 AGeV beam energy. 
$ q depends strongly on the ratio of accepted charged particles to the total number of 

x 40 and 80 AGeV: 7% most central, 158 AGeV: 10% most central 
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charged particles (iV c h(A?/))/(Ar ch ) tot , which is varied by changing the accepted rapidity 
window Ay. However, it is approximately independent of the beam energy. The solid line 
in Fig. |l] labelled $ q c represents the expectation for a system with total net charge zero, 
in which the only correlations are due to global charge conservation. 

$ q c = y/l- (N ch )/(N ch ) tot - 1 (2) 

In order to enlarge any deviations from this trivial effect, the difference A$ q = $ q — $ q c 
is displayed in Fig. It is found that A$ q is close to zero. 

To study the sensitivity of the A<3> q measure a model, describing a quark gluon plasma, 
was investigated 0. This model assumes an ideal gas of massless quarks and gluons in 
equilibrium with zero baryonic chemical potential (/xb = 0). The requirement of entropy 
and local charge conservation during hadronization allows to extract predictions for A$ q 
in different scenarios. Fig. includes two extreme cases: The one labelled "Frozen QGP" 
assumes hadronization only into pions and no diffusion of the net charge in rapidity 
space, so that the initial QGP-like fluctuations are conserved. In fact, this results in 
A$ q values clearly below zero. In the other scenario ("Resonances") hadronization is 
happening entirely into p-mesons. It turns out that in this case the initial fluctuations are 
completely obscured by the subsequent decay of the resonances, which cause a smearing 
of the original QGP fluctuations in rapidity space. 

1.2. (pt) fluctuations 



Figure 3. $ pt for forward 
rapidities (4.0 < y < 5.5) 
as a function of (N pait ). In- 
cluded are p+p, C+C, and 
Si+Si reactions (filled sym- 
bols), as well as centrality se- 
lected Pb+Pb reactions (open 
symbols), all at 158 AGeV. 
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For the study of event-by-event (p t ) fluctuations again the $ x observable, defined in 
Eq. [I], is used. This time, however, x is replaced by the p t of a given particle. Independent 
particle emission will again result in $ pt = 0. The $ pt values, shown in Fig. §, are all 
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corrected for the effects of the two track resolution of the detector. Generally, it is 
found that $ pt is small for all investigated reaction systems: |$ pt | < 10 MeV/c. A weak 
centrality dependence is observable, with a maximum in <3> pt for very peripheral Pb+Pb 
collisions. Also, $ pt is clearly charge dependent: $ pt is always larger for negatively 
charged particles than for positively charged. This may be a reflection of the fact that 
the positively charged particles contain a larger fraction of baryons, which are subject to 
Fermi-Dirac statistics, while the negatively charged particles are dominated by bosons [[J . 
However, $ pt for all charged particles is still higher than for the negatively charged ones, 
indicating that there are additional correlations present. 
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Figure 4. Two particle correlations for p+p, C+C, Si+Si, and central Pb+Pb reactions 
at 158 AGeV. All figures are on the same scale. 

Since $ pt is a global observable, a small value of $ pt does not neccessarily imply the 
absence of any strong correlation. It is also possible that contributions of two effects just 
cancel out each other. Looking at two particle correlations provides a more differential 
way of studying (p t ) fluctuations |7|]. In the procedure employed here, first the p t of a 
given particle is transformed into a cumulative variable x ||: 

Jo* 

x(p t ) = j£ - , where dn/dp[ is the inclusive p t distribution. (3) 
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Then two particle correlation plots are generated by plotting x\ versus X2 for all particle 
pairs inside one event. Figure ^ shows the result for different reactions systems at 158 
AGeV. While there is a clear structure visible in p+p reactions, reflecting the long range 
correlations present in this case, these structures get more and more diluted when going 
to larger systems. On one side this is naturally due to the effect of the increased combina- 
torics between the growing number of particle pairs, an effect that is removed by the $ pt 
measure. On the other side differences in the reaction dynamics between elementary p+p 
and nucleus-nucleus collisions will show up in the two particle correlations. The relation 
between <3> pt and the two particle correlations is still under study. 

Toy model a(T)/T = 1 0% Pb+Pb (5% central) 




Figure 5. Comparison of a temperature fluctuation model (left hand side) with a(T) /T = 
10% to central Pb+Pb data at 158 AGeV (right hand side). Please note that the data 
are on a different scale than in Fig. 

To see how dynamical fluctuations affect the two particle correlations, a model study 
is performed. In this model the only source of fluctuations are the event-by-event fluctu- 
ations of the slope parameter T of the transverse momentum spectra. As can be seen on 
the left panel of Fig. |, T fluctuations of the order of 10% already result in a very promi- 
nent structure in the two particle correlation. A structure on this level is clearly absent 
in the central Pb+Pb data (right panel of Fig. |5|). Here only short range correlations 
(e.g. Bose Einstein correlations) are visible as an enhancement close to the diagonal. A 
comparison of the measured $ pt in the 5% most central Pb+Pb reactions at 158 AGeV 
to a prediction for small T fluctuations || actually suggests that a(T)/T is smaller than 
1%. 

2. Correlations 

2.1. Centrality dependence of h h Bose Einstein correlations 

The radius parameters shown in Fig. ^ are derived from the fit of the Bertsch-Pratt 
parametrization in the LCMS0 

Cbp = 1 + A exp(— i?g ide gg ide — i?^ ut ql ut — Rf ong qf ong — 2Rl ut i ong <?out<?iong) (4) 
2 Longitudinally Co-Moving System. 
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Figure 6. The radius parameters as a function of (N piirt ) for centrality selected Pb+Pb 
reactions at 40 and 158 AGeV. 



to the h~h~ correlation function (c.f.). The Coulomb correction is included in the fit 
procedure, and is applied only to the fraction of real pairs in the c.f. For both beam 
energies the pairs are in the cm. rapidity region 0.0 < y* < 0.5 with (k t ) = 180 MeV/c. 
The radius parameters in side and out direction show a significant increase with the system 
size. .Riong, however, shows no clear evidence for a variation with (A^ part ), except perhaps 
for very central reactions. The results at 40 and 158 AGeV are very similar. 



2.2. Kaon Bose Einstein correlations in central Pb+Pb reactions 
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Figure shows the c.f. of charged kaons together with the applied fit of the Bertsch 
Pratt parametrization. The data points are corrected for the Coulomb interaction and 
the momentum resolution. 

The table summarizes the values of the radius parameters together with statistical and 



systematical errors. Also included are published results from the NA44 collaboration |L0 
on K + correlations that agree quite well with our measurement. 



2.3. Other two particle correlations 

Apart from the effect of quantum statistics, the c.f. reflects also the influence of the 
final state interaction. This can be exploited to gain useful information from correlations 
between non identical particles. Results on 7rp correlations allow to study relative space 
time asymmetries |fTl| , and Ap correlation give access to the source size [0. Additionally, 
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Figure 7. The correlation functions of charged kaons for the 5% most central Pb+Pb 
reactions at 158 AGeV close to mid-rapidity (0.0 < y* < 0.5). 



one can use the c.f. as a tool to learn about the two particle interaction in cases where it 
is unknown, like in the AA case [13]. The above correlations are studied in the variable: 



Q = <JW = 2&* with k* = -{pi — P2) (in pair cm.). (5) 

2.3.1. 7rp correlations 

Since the relative wave function of pairs of non identical particles contains odd terms in 
k* ■ r*, where r* are the emission points in the pair cm., relative space time asymmetries 
should become visible in the ratio R+/R.. Here, R + is the c.f. containing all pairs with 
^out ■ > 0, while R. consists of pairs with k* ut ■ r* < 0, where the out-direction is defined 
in the LCMS. As can be seen from Fig. ^ a clear asymmetry is observed in the data, which 
goes into opposite directions for 7r + p and 7r~p pairs. This mirror symmetry is caused by 
the fact that the aymmetry in R + / R_ is effected mainly by the Coulomb interaction, which 
introduces a dependence on the charge sign. Additionally, Fig. [8] includes the result from 
a RQMD simulation. Due to its long range character the Coulomb interaction is very 
sensitive to the tails in the spatial distribution of the source. The simulated emission 
points have been scaled by a factor of 0.8 in accordance with the analysis of tt + tt~ and 
7r+ / 7r_ P correlation functions. 

2.3.2. Ap correlations 

Figure |] shows the measured c.f., containing 60k pairs with Q < 0.3 GeV/c, together 
with theoretical c.f. fitted to the data. The calculation is based on an effective range 
approximation, using a s-wave scattering length of / = -2.3 fm (singlet) / -1.8 fm (triplet) 
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Figure 8. The ratio R + /R_ 
in out direction for the 20% 
most central Pb+Pb reactions 
at 158 AGeV. 
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The Gaussian source is assumed to be spherically symmetric and static, and is 
defined by the parameter R G . Figure |9| includes two different fits. In the first Rq and the 
A parameter can vary freely. In the second A is fixed to a value that is estimated from 
the background of particle misidentification and the contribution from feed down. Both 
fits suggest a Gaussian source size of Rq = 3-4 fm, which is compatible with the NA49 
result on pp correlations (R G = 4.0 ± 0. 15±g;?l fm) |§. 

2.3.3. AA correlations 

The significance of the measured AA c.f., shown in Fig. 



10, is unfortunately limited by 



low statistics (3500 pairs with Q < 0.3 GeV/c) and does not show any clear structure. 
Nevertheless, it is worthwhile to do a comparison to theoretical expectations, in an at- 
tempt to limit the range of possible parameter values. Therefore, a fit is performed where 
Rg and A are now fixed and the scattering length fo, describing the strengh of the interac- 
tion, is varied. The result indicates that the c.f. would favour a relatively small fo, quite 
independent from the assumed source size. As a comparison a calculation with f = -20 
fm is also shown, which would correspond to the scattering length in the nucleon-nucleon 
case, and which looks rathers unlikely, although is is not totally ruled out. 

2.4. Deuteron coalescence 

From the measurement of proton and deuteron spectra (all at mid-rapidity for central 
reactions^), a deuteron coalescence factor B2 can be derived: 



E, 



d 3 N d 
dp 3 d 



Bo 



P dpi , 



Pd = 2p P 



(6) 



Figure [TTJ displays the extracted B 2 for different beam energies, also including data from 
the AGS and other CERN SPS experiments. B 2 decreases significantly with increasing 
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Figure 9. The Ap correlation function 
for the 20% most central Pb+Pb reactions 
at 158 AGeV. The lines represent fits of 
the calculated c.f. with fixed A parameter 
(dashed) and and free A (solid). 
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Figure 10. The AA correlation function 
for the 20% most central Pb+Pb reactions. 
The lines display the fit results of the cal- 
culated c.f. to the data for different fixed 
Gaussian source radii Rq. 



beam energy. Even in the SPS energy range, taken alone, there is a change by a factor 
of 2. B 2 is also strongly dependent on the system size, as is illustrated in Fig. [12|. Going 
from very peripheral to central reactions, it decreases by a factor of almost 10. This is 
in contrast to the centrality dependence of the radii from Bose Einstein correlations (see 
Fig. |), which can be related to B 2 via fT6|j : 



3 !— n m d 1 



B 2 = -(^fcc) 3 -^ with R G = $Rl dc R long (7) 
4 m n G 



As is demonstrated in Fig. 12 the B 2 derived from the HBT volume changes much less 



with the centrality than the measured deuteron coalescence parameter. 
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Figure 11. The deuteron coalescence factor 
B2 as a function of the beam energy for 
central nucleus-nucleus reactions. 



Figure 12. The deuteron coalescence factor 
B 2 as a function of the system size for min- 
imum bias Pb+Pb reactions at 158 AGeV. 
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